Search for: All records

The process of self-assembly of biomolecules underlies the formation of macromolecular assemblies, biomolecular materials and protein folding, and thereby is critical in many disciplines and related applications. This process typically spans numerous spatiotemporal scales and hence, is well suited for scientific interrogation via coarse-grained (CG) models used in conjunction with a suitable computational approach. This perspective provides a discussion on different coarse-graining approaches which have been used to develop CG models that resolve the process of self-assembly of biomolecules. 

Abstract The formation of biomolecular materials via dynamical interfacial processes, such as self-assembly and fusion, for diverse compositions and external conditions can be efficiently probed using ensemble Molecular Dynamics (MD). However, this approach requires many simulations when investigating a large composition phase space. In addition, there is difficulty in predicting whether each simulation will yield biomolecular materials with the desired properties or outcomes and how long each simulation will run. These difficulties can be overcome by rules-based management systems, including intermittent inspection, variable sampling, and premature termination or extension of the individual MD simulations. Automating such a management system can significantly improve runtime efficiency and reduce the burden of organizing large ensembles of MD simulations. To this end, a computational framework, the Pipelines for Automating Compliance-based Elimination and Extension (PACE²), is proposed for high-throughput ensemble biomolecular materials simulations. The PACE²framework encompasses Candidate pipelines, where each pipeline includes temporally separated simulation and analysis tasks. When a MD simulation is completed, an analysis task is triggered, which evaluates the MD trajectory for compliance. Compliant simulations are extended to the next MD phase with a suitable sample rate to allow additional, detailed analysis. Non-compliant simulations are eliminated, and their computational resources are reallocated or released. The framework is designed to run on local desktop computers and high-performance computing resources. Preliminary scientific results enabled by the use of PACE²framework are presented, which demonstrate its potential and validates its function. In the future, the framework will be extended to address generalized workflows and investigate composition-structure-property relations for other classes of materials. 

Protein mimics such as peptoids form self-assembled nanostructures whose shape and function are governed by the side chain chemistry and secondary structure. Experiments have shown that a peptoid sequence with a helical secondary structure assembles into microspheres that are stable under various conditions. The conformation and organization of the peptoids within the assemblies remains unknown and is elucidated in this study via a hybrid, bottom-up coarse-graining approach. The resultant coarse-grained (CG) model preserves the chemical and structural details that are critical for capturing the secondary structure of the peptoid. The CG model accurately captures the overall conformation and solvation of the peptoids in an aqueous solution. Furthermore, the model resolves the assembly of multiple peptoids into a hemispherical aggregate that is in qualitative agreement with the corresponding results from experiments. The mildly hydrophilic peptoid residues are placed along the curved interface of the aggregate. The composition of the residues on the exterior of the aggregate is determined by two conformations adopted by the peptoid chains. Hence, the CG model simultaneously captures sequence-specific features and the assembly of a large number of peptoids. This multiscale, multiresolution coarse-graining approach could help in predicting the organization and packing of other tunable oligomeric sequences of relevance to biomedicine and electronics. 

Reconstituted photosynthetic proteins which are activated upon exposure to solar energy hold enormous potential for powering future solid state devices and solar cells. The functionality and integration of these proteins into such devices has been successfully enabled by lipid-like peptides. Yet, a fundamental understanding of the organization of these peptides with respect to the photosynthetic proteins and themselves remains unknown and is critical for guiding the design of such light-activated devices. This study investigates the relative organization of one such peptide sequence V 6 K 2 (V: valine and K: lysine) within assemblies. Given the expansive spatiotemporal scales associated with this study, a hybrid coarse-grained (CG) model which captures the structure, conformation and aggregation of the peptide is adopted. The CG model uses a combination of iterative Boltzmann inversion and force matching to provide insight into the relative organization of V 6 K 2 in assemblies. The CG model reproduces the structure of a V 6 K 2 peptide sequence along with its all atom (AA) solvation structure. The relative organization of multiple peptides in an assembly, as captured by CG simulations, is in agreement with corresponding results from AA simulations. Also, a backmapping procedure reintroduces the AA details of the peptides within the aggregates captured by the CG model to demonstrate the relative organization of the peptides. Furthermore, a large number of peptides self-assemble into an elongated micelle in the CG simulation, which is consistent with experimental findings. The coarse-graining procedure is tested for transferability to longer peptide sequences, and hence can be extended to other amphiphilic peptide sequences. 

Fundamental bacterial functions like quorum sensing can be targeted to replace conventional antibiotic therapies. Nanoparticles or vesicles that bind interfacially to charged biomolecules could be used to block quorum sensing pathways in bacteria. Towards this goal, dendronized vesicles (DVs) encompassing polyamidoamine dendron-grafted amphiphiles (PDAs) and dipalmitoyl- sn-glycero -3-phosphocholine lipids are investigated using the molecular dynamics simulation technique in conjunction with an explicit solvent coarse-grained force field. The key physical factors determining the stability of DVs as a function of the dendron generation and relative concentration are identified. The threshold concentration of each dendron generation that yields stable DVs is determined. Dendrons with lower generations rupture the DVs at high relative concentrations due to the electrostatic repulsions between the terminally protonated amines. Whereas, dendrons with intermediate generations demonstrate a mushroom-to-brush transition. Conformational changes in the dendrons expand the outer DV surface, resulting in instability in the DV bilayer. DVs encompassing dendrons with higher generations incur stresses on the bilayer due to their high charge density and spontaneous curvature. The self-organization of PDAs on the DV surface are examined to understand how the asymmetric stresses are minimized across the bilayer. A set of conditions are determined to be conducive for the formation of a single cluster of PDAs that decorates the DV surface like a mesh. Results from this study can potentially guide the design and synthesis of nanoparticles which target quorum sensing pathways in bacteria towards the prevention and treatment of bacterial infections. Furthermore, these nanoparticles can be used in diverse applications in biomedicine, energy or electronics that require synthetic dendronized cells or the adsorption and transport of charged species.